Study of the magnetic anisotropy in a planar model of the La2Cu04 
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We study the anisotropy energy in a planar model for the Cu — O layers of the La 2 Cu0 4 investi- 
gated in previous works. The Tight-Binding (TB) potential of the model was extended out of the 
validity region with the purpose of incorporate it in the potential. Next, the spin-orbit operator 
was considered in the Hartree Fock solution of the full HF problem. It follows that in spite of the 
fact that the extended potential possesses the square symmetry of the crystalline structure, the 
■ anisotropy energy vanishes in its purely 2D formulation. The result indicates that for the prediction 

by the model of the observed non vanishing magnetic anisotropy at zero doping and temperature in 
La2CuOi, its formulation requires of a 3D representation of its Wannier orbitals and the use of d 
like orbitals. The consideration of more realistic 3D character of the model, by including multiple 
CuO planes, also could be important due to the argued absence of long range order in 2D. 
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PACS numbers: 71.10.-w,74.72.-h,74.72.Gh,74.72.Kf,75.10.-b,71.30.+h 



I. INTRODUCTION 



The superconductivity is the intrinsic capacity of certain materials that allows to conduct electric currents without 
resistance and losses of energy in specific conditions. In 1911 it was discovery by Heike Kamerlingh Onnes[l| when 
observed the electric resistance of the mercury disappear abrupt at 4.2K. In the next years several materials were 
discovered showing the occurrence of this phenomenon below certain critical temperature value Tc. A microscopic 
explanation was not arrived until 1960, when Bardeen, Schriefer and Cooper proposed a successful theory, today 
known as the BCS theory. 

With the discovery of the La^CuO^ superconductor in 1986 at Tc = 30AT, was open a new stage , up to now 
devoted to the obtaining and investigation of such kind of high Tc superconductors. In this case, up to nowadays it 
had not been possible to develop a full explanation for high temperature superconductivity. 

The essential characteristic of the first discovered HTc superconductor LaiCuO^, and various others in its class, 
is that it has a crystalline structure containing Cu02 layers, inside a perovskite structure separated by block layers, 
who play a main role as charge reservoirs. Usually, the block layers are insulators and do not have any contribution to 
the low energy electronics states. On another hand in the Cu0 2 layers, the minimum energy of the electronics state 
are around the Fermi level. 

At half filling (hole concentration p=0), the cuprates are antiferromagnetic insulators (AFI), with a Neel temperature 
Tj- \ near 300 K. When the number of holes grow, the antiferromagnetic phase is rapidly destroyed and the cuprates show 

■ a pseudogap phase. If you continue doping with hole, the metallic phase turns up. 
^vq In 1930 Bloch and Wilson [3j developed the band theory, that explain why several materials are metals and another 

are insulators. However the band theory failed in trying to explain the electronic structure of a large number of 
insulators. In this group were the NiO, CoO and La 2 _ x Sr x Cu0 4 , which were predicted as having metallic behavior. 
This result opened the doors for the profound study of the strongly correlated materials, which are characterized by 
showing strong correlation properties among the electrons. Two of the founding fathers of these studies Mott and 
Hubbard, attributed the insulator behavior to the electron-electron correlation. 

Nowadays, a generalized criterion exists: For getting behaviors like the ones of insulators of Mott, it is necessary to 
take into account the short range correlations between electrons of opposite spin which are included in the phenomeno- 
logical Hubbard models [1, @ . It is also clear that the HF approximation, when is considered from first principles in 
a elemental electronic system, should not be able to predict the correlation between opposite spin electrons. In the 
recent works @, 0] by using a one band model, solved in the HF approximation, it was obtained a gap at T = for the 
La 2 Cu0 4 , then predicting its insulator behavior. This reflect that some of the important properties of these materials, 
that are products of the strong electronic correlations, can be predicted by the model in the HF approximation. This 
outcome does not result strange, after taking in mind that, by example the antiferromagnetism, that is a strong 
correlation property, can be derived form a HF solution of the Hubbard model. The main point to note in this, is 
that the considered model is not the original first principle electronic system. In fact the model only differs from a 
Hubbard one in that the near neighbor approximation was not taken. 

Since the Hamiltonian adopted in the model does not includes the spin, the HF solution do not determine a direction 
for the antiferromagnetism. This fact define the basic motivation of this work. It consists in including the spin-orbit 
interaction in the starting Hamiltonian and then study the effect of the spin-orbit interaction in the solutions. The 
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expectation was to study the magnetic anisotropy of the model. Therefore, he specific objectives of our work are the 
following ones: 

1 . Obtain the form of the spin-orbit operator in the frame of the mentioned model for the superconductor material 

2. Afterwards, to evaluate the anisotropy energy of the antiferromagnetic HF states of the model, in the no doping 
limit of the La 2 Cu0 4 . For this purpose the mean value of the spin-orbit operator of the electrons with the 
crystalline field @, 0] is evaluated in the HF states. Since the HF solution is degenerated with respect to spin 
rotations, to evaluate the mean value, the HF states were arbitrarily rotated in the spin structure of their many 
particles 0,0]) by unitary transformations. These spin rotations are the ones representing corresponding space 
rotations around an arbitrary spatial axis. 



II. LA 2 CU0 4 SYSTEM 

In general, the HTSCs presents a tetragonal structure, and have one or more planes Cu0 2 in its structure separated 
by layers of another atoms (Ba, O, La,...). In the Cu0 2 planes, each cooper ion is coupled to four oxygen ions separated 
an approximate distance of 1.9 A. The critical temperature depended of the maximum number of adjoint layers Cu0 2 
(figure Q}. The electronic behavior of the cuprates is very strongly determined by the electronic density in the two 
dimensional layers Cu0 2 . The basic behavior of the Cu0 2 layers, is common for all the cuprate superconductors. 
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FIG. 1: Dependence of Tc of the HTSC with the number of serial planes Cu0 2 . 



A. Crystalline structure and phases diagrams 

The family of monolayer compounds, La 2 _ x Sr x Cu0 4 , where x indicates the magnitude of strontium (Sr) doping, 
crystallize to the tetragonal structure centered in body (bet), usually named T structure, showed in the figure [2] In 
the La 2 _ x Sr x Cu0 4 the planes Cu0 2 are separated by ~ 6.6 A, separated by 2 planes of LaO, that act like charge 
reservoirs, absorbing electrons of the conductor planes under the doping of the material with holes. 

Exist consensus in consider the electronic structure of the La 2 _ x Sr x Cu0 4 , as quasi two dimensional. This low dimen- 
sionality of this electronic system allows to be treated from another viewpoint. For example collective phenomena, 
such as fluctuations in spin density waves (S.D.W) and charge density waves (C.D.W), could be the origin of the 
superconducting mechanism [H, Q . 

In the figure [3J it is illustrated the phase diagrams in the (x,T) plane for the compound La 2 _ x Sr x Cu0 4 , where T 
is the temperature. Even when the absolute antiferromagnetic order and the insulator Mott state persist for a wide 
range of temperatures, it is very sensible at doping. For example, in the T = limit, it is lost at x = 0.04 doping, but 
in the pseudogap (SG) phases and superconductor (SC) ones, exist experimental evidence of a less extended AFM 
order [llj, [H| ■ AH the ideas mentioned before support the possibility that the AFI state (x=0) play a crucial role in 
guaranteeing first: a correlation between pairs of holes (SG), second: a condensation of these holes toward the SC 
state [13 . 
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FIG. 2: Crystalline structure of the La 2 Cu0 4 



In this compound a structural phase transition occurs. At high temperatures the structure is tetragonal, but for low 
temperature the Cu atoms and the six oxygens around them, stray slowly of their positions, forming an orthorhombic 
structure. In the majority of the theoretical studies about this compound, this little distortion is usually ignored. 

The crystalline structure of the LSCO have simple Cu0 2 planes, while the systems Bi 2 Sr 2 CaCu 2 8+y and 
YBa 2 Cu 3 7 (YBCO) has double Cu0 2 planes. The intermediate layers of LSCO do not have the Cu— O chains 
that affects considerably the Cu0 2 planes, in contrast with YBCO, and the complicate structure that present the 
intermediate layers of the Bi 2212 . 



FIG. 3: Phases diagram of the La 2 _ x M x Cu0 4 . AF=Antiferromagnetic, PM=Paramagnetic, I=Insulator, M=Metal, 
SG=SpinGlass, N=Normal, SC=Superconductor, T=thetragonal structure, 0= Orthorhombic structure. 



The band diagram shown in the figure |4] associated to the La 2 Cu0 4 was obtained by means of techniques of 
LAPWfUj]. The last occupied band is half filled, predicting a metallic behavior, and its form suggests a Tight- 
Binding nature for the electron gas that occupied it within the effective environment of the rest of the bands. The 
less bound electrons to the compound La 2 Cu0 4 , are the ones of the CuJ atoms, that do not have closed their last 
shells (3d) , at difference with the situation with the 2 ions. These electrons are the ones that form the last band 
in the corresponding solid. 

The model proposed in the references @, 0] is described in this section. One of the most used pictures to describe 
the movement of the valence electrons in a solid with a periodic potential is the tight-binding model. This model is 
applicable when the overlapping of the atomic orbitals of the neighbors can be considered small and it can be supposed 
that is not necessary to carry out significant corrections to the atomic wave functions, created by the interaction with 
the neighboring orbitals. The tight-binding approximation provides a reconciliation among the seemingly contradictory 
representations of the very localized atomic levels and on another hand the free electron model [Oj , in which the wave 
functions are lineal combination of planes waves. 



La 2 _ I M I Cu0 4 
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III. THE STARTING THEORETICAL MODEL 



FIG. 4: Bands structure for the La 2 Cu0 4 calculated by Horsch y Stephan et.al. 1993. The semifilled band T — X 

predict a tight-binding behavior in the reciprocal plane 



From the previous comments we can suppose that these electrons are strongly correlated to the base cells Cu0 2 , 
preferably toward the corresponding centers of Cu 15]. Then, the lattice will be supposed as the square lattice of Cu 
atoms located at the CuO planes (see figure [5]). 

The electrons not belonging to the last band and the nuclear charge that neutralizes them, play a double role. 
They act in one hand, as an effective polarizable medium that screens the field produced by a point charge external 
to them, and in second place, due to their spatial distribution and magnitude, they guarantee with their action the 
periodic order of the solid, that is modeled with a periodic potential in the lattice of points. 

The Hamiltonian of the model introduced in the references take the form: 



N 2 N N N N N 

i— 1 i—1 i—1 i—1 j^i 2—1 

-2 _* 

where is the kinetic energy term, W 1 (xi) = Wj(xi + R) is the periodic potential that exercise the electrons in 
the partially filled single band, and R runs over the points of the lattice of Cu (see figure [5]), 

Fb(xi) = J2r I ^"y ifee CXPl ' \£ b -y\^ b ' b p. is the potential of the jellium background correspondingly only to 
the surplus of charges which neutralize the electrons that partially filled the band, 

2 1 

V(xi,Xj) — 4 ^ eeo screened Coulomb potential, 

V(£i) = W 7 (xi) + Fb(xi) is the effective potential that feel the electrons that partially filling the last band, 
Uso — 2 5 s VV(xj) • a x p is the spin-orbit interaction term. 
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FIG. 5: Punctual lattice associated to the CuO planes. To obtain the AFM properties of the electronic gas, and 
more general, to release the symmetries restriction, it is separated the lattice in the 2 represented red 
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A. The Hartree-Fock method 



One of the most useful techniques to find approximate solutions to the Schrodinger 's equation is the molecular 
orbitals approximation. In this approach it is considered that the electrons occupy spin-orbital functions of a given 
energy. This technique is based in the variational theorem, which state that the expected value of H for an arbit rary 
wave function is not smaller that the lowest eigenvalue Eq corresponding to the exact solution of the equation [16|]. 
Then, the best wave function is obtained when the energy calculated by the equation @ has a minimum. 



Then, the problem is reduced to find the solution by finding it form in all the allowed many electronic states, 
under the condition of making minimal the energy of the system. In the HF formalism one keeps each electron in 
a spin-orbital, in the field of the nucleus and the other electrons that are in their respective spin-orbitals (fixed), 
subject to the constraint that the spin-orbitals are orthonormal. After this minimization process, the HF equations 
for individual spin-orbitals result: 



fi<f>i = £i<pi, (3) 

where fi is the Fock operator and £j are the eigenvalues that represent the energies of the spin-orbitals (<j>i). 

The Fock operator for an arbitrary electron 'V (equation 0]) can be separated in three terms like it is shown next: 



iV 



/,(l) = ^°(l)+^[J i (l)-^(l)], (4) 
where the Hamiltonian operator H°(l) is defined by: 



H°(l) = i^ + W^xi) + F b (xi). (5) 
2m 

Here and in what follows, the arguments of the wave functions 'V and '2', indicate the collection formed by the 
spacial coordinates and the spin. In absence of interelectronic interactions ([5]) will be the only operator corresponding 
to the movement of a single electron in the field of the nucleus and the other electrons, that do not pertain to the 
partially filled band crossing the Fermi level. 

The direct interaction operator is 



Jj(^a(l) = U dr 2 ^(2) i -£— ^(2)| a (l), (6) 

where J g?T2 represents the integral for the space coordinates and the sum for the spin coordinates. This operator 
correspond to the average potential of an electron in the state <f>j . 
On the other hand the exchange operator is expressed as: 

= {/ ^^(2) ^^ ^(2)1^(1). (7) 

The equation ([3]) is known as the canonical HF equation and the fact that it corresponds to a non linear eigenvalue 
and eigenfunctions problem, implies that to find the solutions becomes necessary an iterative procedures. This fact 
leads to propose a group of test spin-orbitals that allow to obtain a new group of such orbitals, from which a new 
Fock operator is generated. The process is repeated until the energy and the functions converges, such a procedure is 
named the self-consistent field process. 
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B. The approximate basis functions 



In a general way the spin-orbitals (<f>i) are expanded by means of a well-known group of basis functions associated 
with the geometric positions of the nuclei. This basis should satisfy two practical conditions 

1. Should reasonably describe the physics of the problem. 

2. The integrals that appears in the mathematical treatment should be evaluates with relative easiness. 

The spatial components of the spin-orbitals is expanded in terms of the basis functions generally centered in atoms, 
according to ([8]). The expansion which must be infinite in principle, is truncated using a number of functions M, 
being c V { the coefficients determining each atomic orbital contribution Xu m the linear expansion [l7j 



X] c «*> (8) 



Substituting the orbitals expanded according to © in the HF equation [18], a set of matrix equations is obtained 
which can be: 



FC = SCe, (9) 

where F is the Fock matrix for the F^u — (Xn\f\Xv}i S is the superposition matrix among the elements of the base 
S^v = (XnlXv) an d C is the coefficients matrix c^. 

In the beginnings, the developing of HF techniques for systems of many particles, like bands and polyatomic 
complex calculations, was limited by the memory and speed of the machine calculators of that times. It was necessary 
to devise, on the basis of arguments and credible approaches, for example like the hole's Fermi model [l9j . mean field 
approximation, etc, a reduction of the functional space for searching and of the complications of the calculations. An 
useful assumption is to consider that the solution states of the system of equations (J9j) have quantized spin in the same 
direction in all the space, that is to say: 



(f>£ (x) ti-f(s) type a state, 
<f>±{x) ui(s) type f3 state. 



M*, ') = { ffl o 1: (io) 



If the spatial functions 0£ and are identical, the HF calculations is named as restrictive, otherwise unrestricted 
[20( | . However, due to their definitions, both cases are restrictive. 

It can be said that the HF system of equations (JSJ is written in a form being invariant under rotation, because it is 
not considered the spin quantization in one absolute direction. That is, the states to be found are not a type neither 
/3 one. The first in derive the HF equations in this fully unrestricted form was P. A. M. Dirac in the 1930 [2l|. Work 
in this rotation invariant way wide the functional space being investigated. In consequence it allows to extend the HF 
scheme towards the search of ground states showing magnetic character. This possibility is suggested in references 
0, 0] • These works considered a base that allowed to find HF solutions for the spin-orbitals showing spin quantization 
depending of the spatial position. 

The Tight-Binding basis in the one band approximation, that was used in these works consists of 4 sets of basis 
functions: 

^(x, S ) = J|^( S )^exp( i k.RW)^ (x-R W ), (11) 

b z u"' = a z u°*, (12) 
1 x 2 

<A)(x) = - 7 =exp(-— — ), (13) 



where a <C p, r = 1,2 and a z = ±A. 

The equation © projected in the tight binding basis (fTTj) was written in the form: 
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[H k + x (Jk - K k - F k )].B k -' = ej(k) S k .B k ' ! , (14) 

where i?k is Hamiltonian part that corresponds to the kinetic energy and the total effective potential, 
Jk is the direct interaction term, 

Kk is the exchange term, 

.Fk it is the neutralizing bottom potential, 

Sk is the overlapping matrix. 

Also x = InK^tta p an< ^ ^iOO = ^p-s^k) are the dimensionless parameters in those that a is the characteristic 

radius of the functions ipo, p = 3.8 A is the separation between two Cu near neighbors in the absolute red, h is Planck 
constant, m and e are the mass and the charge of the electron respectively. 

In the low overlapping approximation among nearest neighbors (that is to say among different sublattices) the 
functions of the base that correspond to different sublattices with the same spin quantization are the only one, not 
being orthogonal. However, the orthogonalization of different elements of the same sublattice, as well as the unit norm 
of all elements, are valid. The y>o( x — R^) orbitals will be called as Wannier orbitals and represent the probability 
amplitude to hnd an electron in the surroundings of a point of the lattice. 

IV. MAGNETICS PROPERTIES 

Magnetism arises on the sub-atomic level from localized polarization of the electron clouds of certain atoms arising 
from unpaired electrons. This causes the charge on the atom to have a lattice angular momentum. Any flow of charge 
causes additional physical effects on the surroundings, usually referred to as a magnetic effect. In the case of atoms 
the lattice angular momentum of the charge cloud causes a magnetic field perpendicular to the rotation plane of the 
charge. The magnitude of this magnetic or spin moment is dependent on the species of atoms [22j]. How these spin 
moments interact with each other is critical to how different materials are characterized magnetically. When atoms 
are brought in proximity to each other there is a probability of an electron jumping from one atom to another, known 
as the Heisenberg exchange [4]. This interaction probability can indirectly couple the spin moments of the atoms, 
causing the spin moments to align parallel or anti-parallel. In most materials the spin moments are small and aligned 
randomly, leading to paramagnetism as shown in Figure. In some materials however, specifically transition metals 
such as nickel, cobalt, and iron, the spin moments are large, and align in parallel. This causes a lattice spontaneous 
magnetic moment in the material. 

In La 2 _ x Sr x Cu0 4 case, the use of scattering of neutrons techniques |23| . has contributed important information 
about the magnetic structure of these materials. In the figure [5] the AFM profile of the La 2 Cu0 4 in the copper- 
oxygen is shown. The places Cu contains the magnetic moment, whose address rests approximately to 45 grades 
on the link CuO. It happens that the ions 2 and La^ complete their more external layers, while the CuJ lacks 
an electron for it. Then the ion Cu presents an electron without matching up. It is also the less bound one to the 
arrangement. 




FIG. 6: (a) Magnetic structure of the La 2 _ x Sr x Cu0 4 . Schematic representation of the Cu0 2 planes, the cooper and 
the oxygen are represented by means of hollow and padded circles, respectively (b) Magnetic structure of the 

La 2 _ x Sr x Cu0 4 in its planes Cu0 2 
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V. 



VI. ANISOTROPY ENERGY AND MAGNETOCRYSTALLINE ANISOTROPY 

One of the most basic parameters in a magnetic system is the magneto-crystalline anisotropy, that is, the preference 
of the spin moments to align with some fixed crystallographic axes. This is a consequence of the interaction of the 
local environment though its spin-orbit interaction with the electrons. The simplest form of anisotropy is the uniaxial 
one that consists in that the magnetic moment prefers to align throughout a simple axis e, usually called the easy 
axis. 

Considering the operator Xa=i Uso as an perturbation of the initial Hamiltonian [TJ the anisotropy energy will be 
evaluated in this subsection, in the first correction to the HF energy of the system, according to the corresponding 
formula of perturbation theory. 

£ aniSQtropy (n,0) = ]T < (Plipso^li >, (15) 

k,l 

where the > orbitals are obtained of rotating the orbital solution \<fi k ' 1 > around the axis defined by the unitary 
vector n in an angle 9. To rotate the orbital one should interpret the following procedure: Let us suppose that we have 
an orbital that is eigenfunctions of the operator <j z and we make a rotation around the axis defined by the unitary 
vector n in an angle 9. Then, the z axis will become an z'axis, and the orbital will be transformed in an orbital that 
is eigenfunction of ay according to p^|) . 





(16) 


|05 >=U(a,b)\(f> ktl >, 


(17) 


U{a,b) = {^ t £ ), 


(18) 


= cosi^) ~ m * sin (2^ 


(19) 


. As . . ,(j>, 
■n v sax{~) - m x sm(-). 


(20) 



A rotation of the orbital leads to a rotation of the antiferromagnetic structure of the material at T — 0. This can 
be expressed in the following way [27| : 



M R {x) - ^<05(f)|a|05(f)>, (21) 
k,l 

M R (x) = RJ2< MS)\a\4> k j{x) >, (22) 
k,l 

Mr(x) = RM{x). (23) 

where Mr(x) is the magnetization vector rotated around an n axis, in an angle 9 starting from the initial direction 

of the magnetization vector M(x), and R is the tensor that rotates the magnetization vector in the three-dimensional 
space. 

The quantity i5 an i so t r0 p-a(^, &) can be interpreted in the following way: Let us suppose that the system has its 
antiferromagnetic structure in certain direction and that this structure is rotated around the axis n in an angle 9. 
Then -E an isotrop-a(^) 9) is the correction to the HF energy of the system produced by the spin-orbit interaction, when 
the magnetic structure has that direction (see figure [7]). 



VII. ANISOTROPY ENERGY WITH GAUSSIAN WANNIER ORBITALS 



The anisotropy energy was studied in the initial model. That is, considering as the basis functions, the Bloch 
functions constructed in the sublattices by means of defined by Wannier Gaussian orbitals. It was already seen that 
in this two-dimensional model the spin-orbit operator reduces to the following form: 



FIG. 7: Antiferromagnetic structure in the CuO planes rotated around the z axis in an angle 9. 



Uso = ~ 2 2 C 2 ^ Z -i E ^Pv ~ E vPx)- (24) 

In the following developments, it will be considered that g = — 2 mc 2 » wnere m ls the electron mass. 
The expression for the energy to evaluate is given for: 

E(n 7 0)=Y, < ^i\Uso\^: l >, (25) 

k,l 

where B^' l a are the constants for the expansion of the HF solutions of the original problem 
in terms of the basis functions [6|. 

\<&<"\2, a) >= ^u*»(*)£acr) exp(ik- 8<r))<po@ ~ # (r) )- 
The formula (1251) can be wrote as follows: 



E(n, #) = .9 E E E fl #. B *l/a„«r. {ft, 9)g(r, t), (26) 

k,l r,cr z r,a z 

f a ^M0)=u a >U + a z Uu°% 



^ r '*) = jif E e*-^-^) [d 2 x<p (x-Rt)[E x p y - E y p x ]tp (x- R r ). 
In the nearest neighbor approximation, the term g(r, t) take the following form: 

!J d 2 xipo(x)[E x p y - Eyp x ](p (x) if r=t 

En=-i e ' fcx " P Jv™ <Px<po{x + e x np)[E x p y - E y p x ](f (x)+ (27) 
exp4 fc B n P j d 2 f^ (^ + e y np)[E x py - E y p x ]ipo(x) iir^t 

that can be expressed in more explicit form according to: 



A J d 2 xe 2 ~a? [E y x — E x y]e 2 ^ if r=t 

lfc *"P / d 2 fe" ' x+ 2 eH a 7' [^35 „ E x y]e~^> if r^t 



(28) 



where ^n(^) = / n exp(— ^-y) and A 



— Hi 
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The integration regions are taken in what follows of fixed dimensions because the Wannier orbitals are very much 
localized in the Cu lattice centers, (a = 0.25 p, p — 3.8 A ): 
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As the Wannier orbitals are two-dimensional Gaussian functions, they coincide with the wavefunctions of the 
Harmonic two-dimensional oscillator. Therefore, from the analysis already done, we can easily derive the effective 
potential in the vicinity of the Cu as a parabolic potential (figure [5]). 




FIG. 8: Parabolic potential (up) and effective potential (down) that coincide in the vicinity of the Cu atoms. The 

unit of energy is E = 8.3 eV and the length unit p = 3.8 A . 

The components E x , E y of the electric field are derived form the effective potential, see figure [9] They are odd 
functions in the x, y coordinates, a property that makes that the integral factor E y x — E x y is an odd function. As 
the integration regions are symmetrical regarding, at least one of the coordinates, it follows that g{r,t) = 0. 




FIG. 9: Effective potential with the order symmetry equal to 4 of the crystalline lattice, with parabolic form around 
the Cu atoms, Gaussian Wannier orbitals. The unit of energy is E — 8.3 eV and the length unit p = 3.8 A . 

This result evidences that the Gaussian Wannier orbitals in our small overlapping approximation, predict a vanishing 
anisotropy energy. This means that the first correction to the HF energy for the system, determined by the spin-orbit 
interaction doesn't determine an easy magnetization axis. This result is not coinciding with the observations, which 
indicates the existence an easy axis at zero doping and T — 0. However, this difference can be attributed to various 
factors. One of them is the known absence of long range order in the pure 2D systems. Another cause, can be the 
fact that the model slightly simplifies the microscopic structure of the electron system by assuming it as pure two 
dimensional and moreover, by considering Wannier orbitals not resembling open the shell properties of the Cu atoms. 

The reduced form of g(r, t) (|27p indicates that their value only depends on the properties of symmetry of the electric 
field and the orbitals in a vicinity of the Cu, and not in the whole lattice. It points that out is fundamental to more 
deeply formulate the starting basis of states for the model to predict the magnetic properties. 
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VIII. ANISOTROPY ENERGY WITH WANNIER ORBITALS BEING ODD OR EVEN FUNCTIONS 

In this section, it is considered a modification of the orbitals of the model in reference [1, 0] ■ It is supposed that 
the Wannier orbitals are three-dimensional. Their symmetry is assumed with the same properties as the original 
problem. The Wannier obitals in general sense will be denoted by ip(xi, x^, x 3 ). It is supposed that they are even 
or odd functions regarding the inversion of anyone of their coordinates, referred to the Cu atoms in that they are 
centered. In this case the spin-orbit operator should have the usual form in three dimensions : 

Uso = ~ 2~1T2 [<?x-(d y Vp x - dzVpy) + cf y .(d z Vp x - d x Vp z ) + d z (d x Vp y - d y p x )}. (29) 

From the definition of the anisotropy energy (fTSj) and (|2"9"|) . it follows that it is a lineal combination of terms of the 
form: 



E(n, 0) — ^ (AijHijQ + Bij + Hiji + + Bjj_Uyi_ + ( ',, II, ,2 + (',, II, ,2 I, (30) 

where the Hijk± are of the following form: 

HijO = J Vo dxidx2dx 3 ip(xi,X2,X3)E Xz pjip(xi,X2 1 x 3 ) 1 

Hiji± = J Vi ^dxidx2dx 3 ip(xi±p,X2,x 3 )E Xt p : jip(x 1 ,X2±p,x 3 ), 

H lj2 ± = fy 2 ^dx 1 dx 2 dx 3 'ip(xi,X2±p,x 3 )E Xz p : jip(xi±p,X2,x 3 ), 

where the regions of integration V\, V 3x ^ 1 V2y^ can have in general sense three possible spatial symmetries, when are 
defined by: 

V = (x,y,z)\ - § <x,y < f,-| < z < |, 

Vi+ = (x,y,z)\ -§<£< §§,-§ <»<§,-§<«<§, 

V l - = (x,y,z)\-te<x<$,-$<y<z,-§<z<§, 

V 2+ = (x,y,z)\ - 2 <x < |, -| <y < ||,-| < z < | 

' V 2 - = (x,y,z)\ - § < x < §,— ff < y < |,-| < z < |. 

Next it follows that: 



dxidx2dx 3 ip(xi,X2,x 3 )E Xi pjip(xi,X2,x 3 ) =0, (31) 
Vb 



Vi 



dxidx 2 dx 3 ip(xi ±p,X2,x 3 )E Xz p J ip(x 1 ,X2±p,x 3 ) =0, (32) 
dxidx2dx 3 ip(xi,X2±p,x 3 )E Xt pjip(xi±p,x 2 ,x 3 ) =0. (33) 



T 



It can be underlined that the two functions p(x\, X2, x 3 ) and t(x\, X2, x 3 ) have the same parity regarding Xi if both 
are even or odd. 

It is considered in our analysis that the effective potential near the Cu atoms of our model is even. Hence it follows 
that the fields E Xi — —d Xi V are odd functions regarding Xi and once fixed Xi , is even regarding the others coordinates 
x j,k- 

Let us show the (|3HI33p properties now. 
Proof of EU 
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1. The functions i/K^ii x 2, £3) and d Xj tp(xi, X2, X3) has the same parity regarding Xi and using the initial premise 
about the fields, the integrand is an odd function concerning the Xi coordinate. 

2. Let us use the result: If the integral of a function converges absolutely, in an integration region in R 3 that has 
parallelepiped form, then the iterate integrals can commute the integration order (Theorem of Fubini) . 

3. Using the previous result we concluded that J v dxidx2dx3ip(x\, X2, X3)E Xi pjip(xi , X2, X3) = 
Proof of El 

1. If Xi = x\, then tl){x\, X2, X3) and d Xj ij){xi ± p, x%, X3) have different parity regarding Xj and E Xi and continues 
being even regarding Xj, of the initial premise, for what (|32[) is kept. 

2. If Xi — X2,3, 4>(x\, X2,xs) and d Xj ip(xi ±p, X2, X3) have identical parity regarding Xi and using the initial premise 
of the fields, then the integrand is odd regard the Xi coordinate, for that (|3"2")l is kept. 

3. In both cases the before cited Fubini property was employed. 
Proof of M 

1. If Xi — X2, it follows that ij)(xi, X2, X3) and d Xj ip(x\, X2±p, X3) have different parity regard Xj, and E Xi continues 
been even regard Xj, of the initial premise, for what it vanishes. 

2. If Xi = £1,3, tp(xi, X2, X3) and d Xj ip(xi, X2 ±p, X3) have identical parity regard X{ and using the initial premise 
of the fields, then has an odd integrand regard the Xi coordinate, and ([3"3"]l is kept. 

3. In both cases the result of Fubini also was used. 

The results are summarized in the following form: If the Wannier orbitals that we propose are even or odd functions 
regard the space coordinates and the effective potential that it is assumed in the vecinity of the Cu atoms, is an even 
function regards the space coordinates, the anisotropy energy in the model also vanishes and therefore an axis easy 
magnetization is not predicted. Notice that the condition, that Wannier orbitals that we propose are even or odd 
functions, is enough but not necessary condition for the vanishing of the energy. This fits the possibility that using 
orbitals that neither be even or odd, a non vanishing anisotropy energy could follows. 

As it is known that the less bound electrons of the Cu^~ occupy the d orbitals, and the electrons Cu 2 + considered 
into our model, are representing the one partially filling the most energetic band in La,2CuO<i, we can expect that 
the Wannier orbitals of our problem should more closely resemble the <ij,2_ y 2orbitals of Cu. These d x i_ y i orbitals 
are solutions of the hydrogen atom,. Then, we could assume as the potential in the vicinities of the Cu a Coulomb 
potential alternatively as it was considered in the initial model. The Wannier orbitals of the d x 2_ y 2 type are even 
functions regarding the space coordinates. The components of the electric field E Xi are odd functions because derived 
of the coulomb potential. Then we uses the previous result to conclude that, if we suppose orbitals of the d x i_ y i 
type then the model again will predicts anisotropy energy null and non easy axis of magnetization. However, the 
possibility of employing a general linear combination of the d orbitals of Cu in addition to a more realistic 3-D character 
of the model by including multiple CuO planes, opens the the opportunity to describe the observed antiferromagnetic 
anisotropy at low doping and temperatures. 

IX. SPIN-ORBIT OPERATOR IN THE HARTREE-FOCK APPROXIMATION 

To study the influence of the spin-orbit operator in the Hartree-Fock approximation the operator is projected in 
the basis of the orbitals proposed in the model: 

[H k + x (Jk - K k - F k )].B k >< = r,(k) S k .B k '. (34) 
We can include the spin operator to obtain the next matrix equation 



[H k + U k + x (Jk - K k - F k )].B k < ; = ei(k) S k .B k >', 
where the derivation of the term U k in the model is presented next: 



(35) 



13 



<<P^ a '\x,s)\U so \<l>k,i > = <$^\g,8)\M a (EJ y -E v p x )£B!l t \ V £'' , ' ) >, 
< <p£ a *\x,s)\U S o\<t>k,i > = E < ^\x,s)\Aa z (E x p y - E y p x )\^ ] > 
where A= -5^7 

?7 k so - A < <p£< a '\Z, s)\a z .(E x p y - E y p x )\ip^ >, 



U£° = A^a z 5 a ^ E e*^ W - a(r) ) | ^(af- R [t] ){E x p y - E y p x )<p Q (x - flW). 
In the near neighbor approximation the term £7^° take the following form: 



h = j Vt d 2 xipo(x)(E x p y - E y p x )(fo(x), 

It2 = J v 3 d 2 x(p (x =F pe x )(E x p y - E y p x )ip (x), 

7 T3 = Jv , d 2 x(p (xTpe v )(E x p y - E y p x )(p (x). 

For the same reasons that in section (anisotropy energy with orbitals of Wannier gaussian) the five integrals are null. 
h = I T 2 =I T 3 = 0=^ Ug° = 0. 

We can conclude that the spin-orbit operator is null in the matrix Hartree Fock representation, that implies no 
change in the Hartree Fock solution. 

The same analysis could be done with the three dimensional Wannier orbitals and the result using the sec- 
tion (Anisotropy energy with orbitals of Wannier that are odd or even functions) would be the same as for bidimensional 
Gaussian Wannier orbitals. 



X. SUMMARY 

We have studied the magnetic anisotropy properties of La 2 Cu0 4 at zero temperature, in the context of a two- 
dimensional model of the layers CuO, which Hartree Fock solutions were obtained in references @. For this purpose 
the spin-orbit interaction was considered as a perturbation to the initial Hamiltonian of the system and the energy 
was evaluated as a function of the direction of the antiferromagnetic structure of the material, when it is rotated. The 
basic conclusions of the work are: 

1) The magnetic anisotropy energy strongly depends on the space dependence of the orbitals of Wannier and the 
effective potential that the electrons feel in the proximity of Cu. 

2) For the Gaussian Wannier orbitals and their parabolic corresponding effective potential employed in the model, 
it was obtained that the magnetic anisotropy energy vanishes. Then, an easy axis of magnetization doesn't appears 
and the antiferromagnetic structure doesn't align in any direction. 

3) The use of more general Wannier orbitals, as can be Gaussian atomic three-dimensional or orbital d type, 
assuming that they are even or odd, all produce vanishing vanishing anisotropy energy. This follows mainly because 
the orbitals and the potentials are even or odd functions with regard the coordinates. 

4) The approach employed, that only considers non null overlapping among orbitals that are centered in near 
neighbors, is central in the previous conclusions. 
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5) The fact that the lattice of the Cu atoms is squared, strongly influences these results of zero magnetic anisotropy 
energy. 

6) Wannier orbitals that are not spatially even neither odd, can predict non zero values of the magnetic anisotropy 
energy. 
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